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We present our results on the non-perturbative evaluation of the renormalization constant for the 
quark field, Z c/ , in Landau gauge within RI-MOM scheme. Using three lattice spacing we are able 
to isolate lattice artefacts of various origin, both perturbative and non-perturbative. In particular, 
the existence of the dimension-two gluon-condensate is discussed, and confirmed. 
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1. Introduction 

As any other way of computing matrix elements, lattice QCD requires calculation of the renor- 
malisation constants, which relate bare quantities to the physical ones. In our particular case the 
renormalization constants will be dependent not only on the scale, but also on the lattice param- 
eters, such as light quark mass and lattice spacing. This dependence turns out to be remarkably 
strong in case of the latter. Even if the lattice computation contains only 0(a) lattice artefacts, the 
bare quantities differ from the continuum ones by 0(g 2 ) ~ 0(l/log(a 2 )) which is unacceptable. 
Renormalization restores the 0(a 2 ) accuracy. While some perturbative methods to calculate the 
constants exist, it is obviously preferable to do it in a consistent manner, i.e. non-perturbatively. 
Although in some cases the results may be similar, there is no guarantee for the validity of the 
perturbative calculation. 

Our method of preference is the so-called RFMOM scheme. It involves the computation of 
Green functions of quarks, gluons, ghosts, at large enough momenta in a fixed gauge, usually the 
Landau gauge. This gives the renormalisation constant Z(ju) at many values of the scale pL. We 
have to note here, that even though we advocate the completely non-perturbative approach, the 
conversion to other renormalization schemes has to fall back onto perturbative methods. Here we 
use perturbative QCD to convert MOM into MS and run to 2 GeV. The running of Zmom(m) is a 
very powerful testing tool: perturbative QCD is only useful if we are in the perturbative regime, 
i.e. at large enough momenta. The only way to check if this is the case is to compare lattice data 
with the perturbative running. It turns out that this is not always so. 

To calculate the renormalization constants we must fix the gauge, and we perform all calcula- 
tion in the Landau gauge. Wilson operator expansion suggests the presence of the non-vanishing 
vacuum expectation value of the only dimension-two operator in the Landau gauge: A 2 = A^A"^, 
and that it is not small [1, 2], [3]. This contribution to the OPE will scale as l/p 2 up to logarithmic 
corrections. To argue that this is a continuum feature we must check that the 1 / p 2 term scales well 
with the lattice spacing when expressed in physical units. 

In general, our museum of artefacts has a number of exhibitions. The 0(a 2 ) artefacts can be 
quite large since we consider large momenta, while finite volume artefacts are minor. There are 
two main types of 0(a) artefacts: 0(a 2 p 2 ) artefacts which respect the continuum 0(4) rotation 
symmetry, and the ones which do not. The latter appear due to the explicit breaking of the rotational 
symmetry, which, on the lattice, is reduced to the hypercubic symmetry H4. After removing these, 
we will identify the 0(a 2 p 2 ) artefacts non-perturbatively by doing a fit of the running Z(ju) which 
will include the perturbative running, the condensate and a term proportional to a 2 p 2 . 

For the elimination of hypercubic artefacts several methods have been proposed in literature: 
the democratic one, the perturbative correction and the non-perturbative "egalitarian" one. The 
perturbative one deserves a separate paper, so we will focus on two others. We will demonstrate to 
what extent the "half-fishbone" structure, which raw lattice results for Z q always exhibit and which 
is a dramatic illustration of hypercubic artefacts, is corrected by every method. 

Although the issues raised here concern all the renormalisation constants as well as the QCD 
coupling constant, we will concentrate in the following on Z q , that renormalizes the quark field 

q R =Z l q ,2 q B (1.1) 
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where qg (qn) is the bare (renormalized) quark field. 

2. Simulation details 

For the details about twisted mass and tree-level improved Symanzik gauge actions please 
refer to refs. [6, 7, 8, 9]. Here we will just briefly mention the essential parts for our presentation 
(See Tab. 1). 

The Wilson twisted mass fermionic lattice action for two flavours of mass degenerate quarks 
reads (in the so called twisted basis [4, 10] ) 

Sfm = a4 £ { Xx [Dw + mo + ijs Wq] Xx } , 

x , (2.1) 

(XT 

£>w = 27m(V„+V*)- y V m V*, 

where mo is the bare untwisted quark mass and \i q the bare twisted quark mass, T3 is the third Pauli 
matrix acting in flavour space and r is the Wilson parameter, which, as usually, is set to r = 1 in the 
simulations. The twisted Dirac operator is defined as 

Aw = D W + mo + z'y 5 T 3 (2.2) 

The bare quark mass mo is related as usual to the so-called hopping parameter K, by K = 
1/(8 + 2amo). Twisted mass fermions are said to be at maximal twist if the bare untwisted mass is 
tuned to its critical value, m cr i t . This corresponds to setting the so-called untwisted PCAC mass to 
zero, so it requires additional tuning during the simulation. 

In the gauge sector the tree-level Symanzik improved gauge action (tlSym) [5] is used: 

^ = fE(*0 I {l-ffieTr^^l+fcj £ {l-MeTr^)}), (2.3) 

where j8 = 6/gg, go being the bare lattice coupling and we set b\ = —1/12 (with bo = 1 — 8Z?i as 
dictated by the requirement of continuum limit normalization). The overview of the ensembles used 
for this work can be found in Tab. 1. From our previous simulations we deduced that fortunately 
Z q does not depend within uncertainties on the sea-quarks mass and on the lattice volume. So here 
we keep both quark mass and volume fixed. 

To compute the renormalization constants for the quark propagator we need to fix the gauge 
and calculate the 2-point quark Green functions on each configuration. We do this using a local 
source taken at a random point on the lattice which reduces the correlation between successive 
configurations: 

S{y,xo)t a ' MA =D^{y,x) a ^so b ^{x) so\*{x) = 8 Vo 8 b , bo 8 pA (2.4) 

where i = u,d labels isospin. Thereafter we perform the Fourier transform of the incoming quark 
which is a 12 x 12 complex matrix 

SiW^e-^y-^St&xo) (2.5) 

y 
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p 


a\l q 


Volume 


# confs dim 3 


# confs dim 4. 




0.004 




120 


100 


3.9 


0.0064 


24 3 x 48 


20 


20 




0.010 




20 


20 




0.006 


24 3 x 48 


100 


100 




0.003 


32 3 x 64 


20 




4.05 


0.006 


32 3 x 64 


20 






0.008 


32 3 x 64 


20 






0.012 


32 3 x 64 


20 




4.2 


0.0065 
0.002 


32 3 x 64 
24 3 x 48 


100 


100 
100 



Table 1: ETMC ensembles used for this work, "dim 3" refers to the local operators (V, A, S, P) while "dim 
4" refer to the operators with one derivative (044, 014, etc). 

and the outgoing quark 

s] 5 ( P ) = y 5 st(p)y 5 (2.6) 

where u = d;d = u. We define the quark renormalization constant Z q as 

Z M = ^ <TrlS- 1 (p)p/]> (2.7) 
where < ... > means here the average over the chosen ensemble of thermalized configurations. 

3. Artefacts 

3.1 Hypercubic H (4) -extrapolation 

A first kind of artefacts that can be systematically cured [11, 12] are those arising from the 
breaking of the rotational symmetry of the Euclidean space-time when using an hypercubic lattice, 
where this symmetry is restricted to the discrete H(4) isometry group. It is convenient to compute 
first the average of any dimensionless lattice quantity Q{ap tl ) over every orbit of the group H(4). 
In general several orbits of H(4) correspond to one value of p 2 . Defining the H(4) invariants 

P [A] =tpi p [6] =ipt (3-D 

fi=l n=l 

it happens that the orbits of H{4) are labeled by the set p 2 ,a 2 pW,a 4 /?^. Group theory tells us that 
any H (4) -invariant polynome will only depend on the four invariants p^ with/= 1,2,3,4 [11, 12]. 
As a consequence of the upper cut for momenta, the first three of these invariants suffice to label 
all the orbits we deal with and hence any presumed dependence on p^ is neglected. We have 
to admit that our action is obviously non-polynomial, but perturbative results suggest that up to 
a certain order the renormalization constant will be such, so it is sensible to assume this feature. 
Later we will be able to check that it leads to consistent results, proving our assumption. Moreover, 
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in the continuum limit the effect of a 2 p^ 4 \a 4 p^ vanishes. We can thus define the quantity Q(app) 
averaged over H{4) as 

G(oV,«VVp [ Va£cd). (3.2) 

If the lattice spacing is small enough such that e = a 2 p^/p 2 <C 1, the dimensionless lattice 
correlation function defined in Eq. (3.2) can be expanded in powers of e and truncated: 



G(aV A [ V/> [ VaJcd) = Q(a 2 p 2 ,a 2 A 2 QCD ) + ^ 



2P [4] 

a 2 ^ T + --- (3.3) 
£ =o P 



The basic method is to fit from the whole set of orbits, sharing the same p 2 , the coefficient dQ/de 
and get the extrapolated value of Q, free from H{4) artefacts. The fit is performed in a momentum 
window between (p — 8,p + 8) to extract the extrapolated value of Q for the momenta p in the 
window, and then we shift to the next window etc. 

Such method, which we call SWF, "sliding window fit", is quite reliable since the extrapolation 
does not rely on any particular assumption for the coefficients of the Taylor expansion. We can then 
estimate the systematic error by varying the width of the fitting window. If we further assume that 
the coefficient 



R(a 2 p 2 ,a 2 A 2 QCD 



de 

has a smooth dependence on a 2 p 2 over a given momentum window, we can expand R as 



dQ[a 2 p 2 ,0,0,a 2 A 2 QCDj 

e=0 



(3.4) 



R = CaipA + c a \ p A,a 2 p 2 (3.5) 
and make a global fit over a wide range of momenta: 

Q(a 2 p 2 ,a 4 pW , a 6 p® , a 2 A\ CD ) = Q(a 2 p 2 , a 2 A 2 Q CD ) + c a2 ^a 2 ^ + c a 4 p4 a 4 p [4] + • • • (3.6) 
We will refer to this method as OWF, "one window fit". 



4. Lattice results and Hypercubic corrections 

The hypercubic artefacts are clearly visible on the raw lattice data for Z^ att as the so-called 
"half-fishbone" structure [13] shown in fig 1. The color code shows the value of the ratio pW / (p 2 ) 2 
which is between 0.25 and 1 . The values which are closer to one are "least democratic" or "tyran- 
nic" ones. We see, as expected, that the tyrannic points are more affected by the artefact. We also 
see that the gap between Z^ att (a 2 p 2 , a A p^ , a 6 p^ , ap 4 , a 2 Ag CD ) at a given q 2 can be as large as 0.07, 
i.e. about 10%. Taking a naive average without a correct treatment of this artefact would leave a 
systematic upward shift of about 5 %. 

The oldest method to alleviate this problem is the "democratic selection". It amounts to keep- 
ing only, say, the violet points in Fig. 1, or, in case of low statistics, being less restrictive, the violet 
and yellow ones. It works to a certain extent, but obviously we still have an upward shift compared 
to H4-method, which is not negligible. 
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Figure 1: On the l.h.s. we show the raw data for j8 = 3.9, Zf*-(a 2 p 2 ,a 4 p^,a 6 p^,ap4,a 2 A^ CD ). The "half- 
fishbone structure" due to hypercubic artefacts is clearly seen. There is one point for every cubic (3-D) orbit. 
On the r.h.s we show, using the OWF, the non-perturbatively corrected data. As on l.h.s., there is one point 
for every cubic (3-D) orbit. 

While performing systematic treatment using H4-method, we only expand up to pM since the 
higher order terms turn out to be negligible. In the l.h.s of Fig. 1, we show the comparison of 
hypercubic corrected data after applying OWF and SWF for the case j8 = 3.9 . The difference does 
not appear to be large which is rather encouraging. The OWF gives a slightly smoother result. 



J3 


a 1 fm 


C a 2pA 


CaApA 


Ca2pA/g 2 


CaApA/g 2 


X 2 /d.o.f 


3.9 


0.00689 


0.067(4) 


-0.0149(10) 


0.044(3) 


-0.0097(7) 


4.1 


4.05 


0.00456 


0.065(3) 


-0.0144(5) 


0.044(2) 


-0.0097(3) 


0.53 


4.2 


0.00303 


0.055(11) 


-0.0124(4) 


0.039(8) 


-0.0089(3) 


0.98 



Table 2: Results for the slope in a- p W/p 2 andaV 4 ! and the same divided by g 2 in the one window fits. 

The fitted values of c a i P A and c a ApA from the one window fit are given in Tab. 2 as well as the 
same divided by g 2 , since perturbation theory expects at least for c a 2 P A to be oc g 2 _ Before dividing 
by g 2 a small scaling violation is apparent, and much less so afterwards. The % 2 in Tab. 2 is not 
good for j3 =3.9, apparently due to some structure at the lower end of the plot, but remember it 
uses only two hypercubic parameters. Of course, this point is further away from the continuum 
than both other points. 

5. Results and Conclusion 

Once the hypercubic corrections are performed we make a fit for every j8 according to the 
following formula 

yNPsubtracted/ 2 2\ _ 7 , 2 ,,2\ / i , c \overp2 \ 7 2 ( c n 
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where coz f/ (^ 2 ,M 2 ) is the perturbative running of Z q [14], c c a P 2 is a hypercubic insensitive lattice 
artefact, and c\ overp 2 is due to the g 2 (A 2 ) condensate. Combining several analysis methods we get 

g 2 (^l 2 )(A 2 )^ = 2.25(40) (^ 5 ) GeV 2 /I = lOGeV 
ZP ert ((10GeV) 2 ,^ are ) = 0.741(3)(7)-0.305(25)(g 2 are -i-5) 
ZP ert ((2GeV) 2 ,g 2 are ) = 0.773(3)(7)-0.319(20)(g 2 are -l-5) 

(5.2) 

To conclude, we calculated Z q from the ETMC gauge configurations in the RI-MOM scheme. 
We demonstrated effectiveness of the H4-method of removing hypercubic artefacts. From the 
resulting hypercubic corrected function Z q (a 2 p 2 ,a 2 Aq CD ) we perform a fit to an Ansatz consisting 
of the perturbative running, a non perturbative 1 / p 2 term and rotationally-symmetric lattice spacing 
artefact proportional to a 2 p 2 . The fits are good and a p term scales almost perfectly in lattice 
units, as expected. The g (A 2 ) term scales rather well in physical units as expected, though less 
accurate. Therefore we advocate the presence of the non-perturbative condensate in the quark 
renormalization constant. More detailed analysis and even more rigorous proof will be presented 
in a separate paper soon. 
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